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Abstract 

We analyze the results for infinite nuclear and neutron matter using the standard relativistic 
mean field model and its recent effective field theory motivated generalization. For the first time, 
we show quantitatively that the inclusion in the effective theory of vector meson self-interactions 
and scalar-vector cross-interactions explains naturally the recent experimental observations of the 
softness of the nuclear equation of state, without losing the advantages of the standard relativistic 
model for finite nuclei. 
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I. INTRODUCTION 



In the quest for a unified model describing both finite nuclei and nuclear matter properties, 
the relativistic mean field (RMF) approach to quantum hadrodynamics has become a very 
popular tool. The original linear a-u model of Walecka was complemented on an empirical 
basis with cubic and quartic non-linearities P] of the a meson to have proper results. We 
shall refer to this approach as standard RMF model. In recent years novel approaches based 
on modern concepts of effective field theory (EFT) and density functional theory (DFT) 
for hadrons, in different implementations, have been developed for the relativistic nuclear 
many - b0dy p roblem B fl S Q fl 3. One of the g oa.s is to overeon. the a d h oc nature 
of the previous models and the presence of uncontrolled approximations. Here we address 
the chiral effective Lagrangian model proposed by Furnstahl, Serot and Tang |3], whose 
mean field treatment is hereafter called E-RMF model. A nuclear EFT contains all the non- 
renormalizable couplings consistent with the underlying symmetries of QCD. The effective 
Lagrangian is expanded in powers of the fields and their derivatives, and the possible terms 
and their importance is systematically categorized P, U, 0, |(| . None of the couplings present 
at a given order can be arbitrarily dropped out without a symmetry argument. 

The free parameters of the E-RMF Lagrangian have been optimized by fitting to the 
ground-state properties of a few doubly- magic nuclei Q]. The truncation of the effective 
Lagrangian at the first lower orders is validated by the fact that when one includes up to 
fourth order terms in the expansion, the E-RMF fits (parameter sets Gl and G2 determined 
in Ref. 3| ) display naturalness and the results are not dominated by the last terms retained 
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121. 



Very recently, the flow of matter in heavy-ion collisions has been analyzed to determine 
experimentally the pressures attained at densities ranging from two to five times the sat- 
uration density of nuclear matter Q]. Nuclear collisions in man-made laboratories can 
nowadays compress nuclear matter to such extreme densities, which occur in Nature within 
neutron stars and core-collapse supernovae. Nuclear models of various sorts that are tuned 
to the known properties of terrestrial nuclei at normal densities or to nucleon-nucleon (NN) 
scattering data are commonly extrapolated for predictions of those very dense systems. The 
recent breakthrough achieved in the experimental determination of the equation of state 
(EOS) of high-density matter ^| motivates us to study the applicability of the RMF and 
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E-RMF models to those conditions. The performance of the E-RMF a ppr oach for finite 



nuclei has already been demonstrated in a series of previous works 

II. DISCUSSION 

The energy density functional of the E-RMF model is written as 
S(t) = 



yyyy. 
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The energy density contains tensor couplings and scalar-vector and vector-vector meson 
interactions in addition to the standard scalar self-interactions K3 and k&. One interprets 
the E-RMF formalism as a covariant formulation of DFT. That is, the mean field model 
approximates the exact (but unknown) energy functional of the many-nucleon system, which 
includes all higher-order correlations, by expanding it in powers of the meson fields. The 
latter play the role of auxiliary Kohn-Sham potentials. The unresolved dynamics beyond 
mean field is encoded in the fitted coupling constants of the effective theory and the intro- 
duction of new interaction terms aims at an improved representation of the nuclear energy 
unctional. Further insight into the concepts of the E-RMF model can be gained from Refs. 
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For uniform nuclear matter all of the terms with gradients in the energy density £ vanish 



and the nucleon density is given by p = -^kp (7=4 for symmetric nuclear matter, 7=2 
for neutron matter, and hp is the Fermi momentum). The pressure P follows from the 
derivative of the energy density with respect to t ie nucleon density: P = p 2 J^(£/p). 

The recent experimental observations y, ^| rule out any strongly repulsive nuclear 
EOS. The constraints on the EOS of symmetric nuclear matter at zero temperature, derived 
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FIG. 1: (a) Zero temperature EOS for symmetric nuclear matter. The shaded area represents the 
region consistent with the experimental data The nuclear matter incompressibility of each 

interaction is expressed in MeV. (b) Zero temperature EOS for neutron matter. The upper and 
lower shaded areas depict the regions compatible with experiment after inclusion of the pressure 
from asymmetry terms with, respectively, strong and weak density dependences ["L^ . 

experimentally by analyzing elliptic and transverse flow observables in nuclear collisions, are 
shown in Fig. la along with the predictions from different RMF and E-RMF sets. We can 
see that the calculations of dense matter based on NL1 and NL3 [JjJ deviate drastically 
from experiment, while the E-RMF calculations with Gl (to some extent) and G2 (to a 
better extent) agree more with the allowed region. A similar situation prevails in the EOS 
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of neutron matter (Fig. lb). The figures show that conventional sets such as NL1 and 
NL3 are not well suited for describing the EOS of dense matter, principally due to the 
fact that without the additional self- and cross-interactions the density dependence of the 
mesonic mean fields is too poor Q]. Notice that although the incompressibility of 
NL1 is = 212 MeV, which is within the empirical limit, the resulting EOS is stiff at 
higher densities and does not follow the experimental trend with increasing density. On the 
other hand, the E-RMF calculations explain better the situation without any forced changes 
in the fitted parameters or the formalism, and the consistency of G2 with experiment is 
outstanding. As depicted in Fig. 1, the microscopic Dirac-Brueckner-Hartree-Fock (DBHF) 
theory jl^, which starts from the NN interaction in free space, also suggests an EOS for 
dense matter much softer than in the usual RMF approach. One notices that the DBHF 
results, for both nuclear and neutron matter, lie within the regions of the EOS compatible 
with the boundaries extracted from experiment. 

The success of G2 in nuclear matter is possible owing to the bulk couplings Co, Vii an d V2 
which confer an extra density dependence to the scalar and vector self-energies. If using the 
functional form (pQ) one attemps to get a soft EOS for dense matter with vanishing rji and 
772, for realistic saturation conditions, this tends to enforce large unnatural values of Co- If 
with 771 = 772 = one keeps (0 at acceptable values, then K4 tends to become negative. With 
inclusion of the additional parameters (0, rji and r] 2 one can agree better with experiment, 
have K4 > and a not very large (q value, as in the case of G2. The experimental data of 
Fig. 1 also are explained reasonably well Q] by the calculations of Akmal et al j^J which 
employ the Argonne v is interaction, though the use of this interaction for finite nuclei is still 
limited. It will be worth exploring the situation for relativistic models that taking advantage 
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of density-dependent coupling vertices are consistent with the DBHF theory 

In Table I we present the results for the mass and radius of a neutron star. The structure 
of the neutron star is determined through the Tolman-Oppenheimer-Volkov equation, which 
follows from hydrostatic equilibrium in strong gravitational fields and general relativity 22] , 
combined with the EOS for -P(p). One sees that the neutron star properties computed 
with G2 are very much in accord with the available experimental information. The critical 
temperatures for the liquid-gas phase transition from our calculations of nuclear matter 
at finite temperature are also displayed in Table I. These values are governed by the EOS 
around saturation where the E-RMF sets produce similar results to NL3. Hence, one expects 
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TABLE I: The radius R (in km) and the mass ratio M/Mq of a neutron star, and the critical 
temperature T c (in MeV) for the liquid-gas phase transition in nuclear matter and in asymmetric 
matter with asymmetry a = (p n — p p )/(p n + Pp) = 0.2. 
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12.9 


13.7 


13.8 


13.7 





that in the finite nuclei regime Gl and G2 yield results on a par with the celebrated NL3 
parametrization. That this is indeed the case has been proven in previous works and therefore 
we shall not attempt to present more results in this respect here. We rather refer the reader 
to the literature 3, 0, 10, 11, 12, where it has been shown that the bulk and single- 



particle properties of finite nuclei are well described by the E-RMF model with a similar 
quality to the standard RMF parametrizations, not only for stable isotopes but also for 
exotic isospin-rich nuclei far from the valley of (3 stability. 

In choosing between the parameter sets Gl and G2 for further calculations we prefer G2. 
It is worth noting that G2 presents a positive value of the $ coupling constant (^4), as 
opposed to Gl and to most of the successful RMF parametrizations, such as NL3. Though 
the energy spectrum strictly has no lower bound with a negative K4 |23||, such negative value 
is necessary in the standard RMF model to get the results closer to experiment. On the 
other hand, to have a positive K4 it is not indispensable to make two parameter sets (one 
for light and one for heavy nuclei) as in Ref. 



III. SUMMARY 



In conclusion, we have shown for the first time that the E-RMF parameter sets give a soft 
EOS both around saturation and at high densities which is consistent with the measurements 



of kaon production [15( and the flow of matter 14] in energetic heavy-ion collisions, and 
with the observed neutron star masses and radii. With the systematic inclusion of new 
interaction terms under the guidance of EFT techniques, and without forcing any change of 
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the parameters initially determined from a few magic nuclei, the E-RMF calculations with 
the G2 set explain finite nuclei and nuclear matter in a unified way with a commendable level 
of accuracy in both the cases. The authors are not aware that it has been possible to prove 
this for other effective nuclear models. At present, the E-RMF approach can be considered 
as a salient step towards a unified theory for finite nuclei as well as for infinite nuclear 
matter and highlights the potential of the application of effective field theory formulations 
for nuclear structure studies. 
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